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Abstract 



The production of the prompt charm mesons D , D+, D*+, and their antip articles, was measured with 
the ALICE detector in Pb-Pb collisions at the LHC, at a centre-of-mass energy ^/s^ = 2.76 TeV 
per nucleon-nucleon collision. The pt-differential production yields in the range 2 < pt< 16GeV/c 
at central rapidity, \y\ < 0.5, were used to calculate the nuclear modification factor /?aa with respect 
to a proton-proton reference obtained from the cross section measured at y/s = 7 TeV and scaled 
to = 2.76 TeV. For the three meson species, Raa shows a suppression by a factor 3^, for 
transverse momenta larger than 5 GeV/c in the 20% most central collisions. The suppression is 
reduced for peripheral colhsions. 
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1 Introduction 

A high-density colour-deconfined state of strongly-interacting matter is expected to be formed in high- 
energy collisions of heavy nuclei. According to calculations of Quantum Chromodynamics (QCD) on the 
lattice, under the conditions of high energy density and temperature reached in these collisions, a phase 
transition to a Quark-Gluon Plasma (QGP) occurs. In such conditions, the confinement of quarks and 
gluons into hadrons vanishes, and chiral symmetry is restored (see e.g. ||T]-l3|). Heavy-flavour hadrons, 
containing charm and beauty, ai^e effective probes of the conditions of the medium formed in nucleus- 
nucleus collisions at high energy. Hard partons, including gluons, light-flavour quaifc, and heavy quaite, 
are produced at the initial stage of the collision in high-virtuality scattering processes. They interact 
with the medium, and are expected to be sensitive to its energy density, through the mechanism of 
parton energy loss. This QCD energy loss is expected to occur via both inelastic (medium-induced 
gluon radiation, or radiative energy loss) HIS] and elastic (collisional energy loss) processes. In 
QCD, quai^ks have a smaller colour coupling factor with respect to gluons, so that the energy loss for 
quarks is expected to be smaller than for gluons. In addition, the 'dead-cone effect' should reduce 
small-angle gluon radiation for heavy quarks with moderate energy-over-mass values irTUTTI. thus further 
attenuating the effect of the medium. Instead, other mechanisms, such as in-medium hadron formation 
and dissociation |[T2l[T3l . would determine a stronger effect on heavy-flavour hadrons, characterized 
by smaller formation times than light-flavour hadrons. Finally, low-momentum heavy quarks may be 
to some extent thermalized in the hot and dense system through rescatterings and in-medium resonant 
interactions |[T4l . 

One of the observables that are sensitive to the interaction of hard partons with the medium is the nu- 
clear modification factor Raa- This quantity is defined as the ratio of particle production measured in 
nucleus-nucleus (AA) to that expected from the proton-proton (pp) spectrum scaled by the average num- 
ber (A'^cou) of binary nucleon-nucleon collisions occurring in the nucleus-nucleus collision. Using the 
nuclear overlap function, which is defined as the convolution of the nuclear density profiles of the col- 
liding ions in the Glauber model \[151 . the nuclear modification factor of the transverse momentum (pt) 
distribution can be expressed as: 

1 dyVAA/dpt 
(rAA) dapp/d;?t 

where the AA spectrum corresponds to a given collision-centrality class and {Taa) is the average nuclear 
overlap function for that centrality class and is proportional to (A'^coii)- In-medium energy loss detemiines 
a suppression, Raa < 1, of hadrons at moderate-to-high transverse momentum (pt^2 GeV/c). Given 
the aforementioned properties of parton energy loss, in the range pt^lO GeV/c where the heavy-quark 
masses are not negligible with respect to their momenta, an increase of the Raa value (i.e. a smaller 
suppression) is expected when going from the mostly gluon-originated light-flavour hadrons (e.g. pions) 
to D and B mesons (see e.g. II101I16II ): R'^ < R^^ < R^a- The measurement and comparison of these 
different medium probes should provide a unique test of the colour-charge and mass dependence of 
parton energy loss. 

Experiments at the Relativistic Heavy Ion Collider (RHIC) measured a strong suppression, by a factor 
4-5 at pi > 5 GeV/c, for light-flavour hadrons in central Au-Au collisions at ^/s^ = 200 GeV ifTTl . 
An even stronger suppression — up to a factor 7 at /^t ~ 6-8 GeV/c — was observed in central Pb-Pb 
collisions at y/s^ = 2.76 TeV at the Large Hadron Collider (LHC) |[T8ti20l . At RHIC, the suppression 
of heavy-flavour hadrons, measured indirectly from their inclusive decay electrons 112 111221 . was found to 
be compatible with that of pions and generally stronger than most expectations based on radiative energy 
loss II23U241 . At the LHC, a measurement by the CMS Collaboration indicates a strong suppression, by a 
factor about 3, in the nuclear modification factor of non-prompt J/i/a particles from B meson decays 



We present the first measurement of the nuclear modification factor for D*', D+, D*+ mesons, and their 
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antiparticles, in Pb-Pb collisions at ^/s^ = 2.76 TeV, carried out using the ALICE detector. The ex- 
perimental apparatus ll26l is briefly presented in Section |2] where the Pb-Pb data sample used for this 
analysis is also described. The D meson signals are extracted using a selection based on displaced decay 
vertex reconstruction and particle identification of the decay products, as presented in Section |3] The 
corrections applied to obtain the pt-differential production yields, and the estimation of the systematic 
uncertainties are described in Sections |4] and |6] respectively. The production of D mesons was measured 
in proton-proton collisions at ^/s = 7 TeV and compared to perturbative QCD (pQCD) predictions ll27l . 
The reference for the /?aa measurements was obtained by scaling these results to the Pb-Pb energy via 
a pQCD-driven approach and was validated by comparing to data from a limited-statistics pp sample at 
this energy ll28l . This is discussed in Section |5] The results on the D°, D+, and D*+ nuclear modifi- 
cation factors as a function of transverse momentum and collision centrality are presented in Section |7] 
The results are compared to the charged particle Raa measured with the ALICE detector (19], to the 
non-prompt J/ i/a results by the CMS Collaboration ||25| . and to model predictions. 

2 Experimental apparatus, data sample, event reconstruction and selection 

The ALICE detector, described in detail in Il26l . consists of a central ban^el composed of various detectors 
for particle reconstruction at midrapidity, a forward muon spectrometer, and a set of forward detectors for 
triggering and event characterization. In the following, the subsystems that are utilized in the D meson 
analysis will be briefly described. In particular, the Inner Tracking System (ITS), the Time Projection 
Chamber (TPC), and the Time Of Flight (TOE) detector provide charged particle reconstmction and 
identification in the central pseudo-rapidity region < 0.9). They are embedded in a 0.5 T magnetic 
field parallel to the LHC beam direction (z-axis in the ALICE reference frame). The VZERO detector 
and the Zero Degree Calorimeters (ZDC) are used for triggering and event selection, and the TO detector 
to measure the start time (event time-zero) of the collision. 



The data from Pb-Pb collisions at centre-of-mass energy ^/s^ = 2.76 TeV used for this analysis were 
recorded in November and December 2010 during the first run with heavy-ions at the LHC. The events 
were collected with an interaction trigger based on the infomiation of the Silicon Pixel Detector (SPD) 
and the VZERO detector. The SPD is the innennost part of the ITS. It consists of two cylindrical layers 
of siUcon pixel detectors located at radial positions of 3.9 and 7.6 cm from the beam line, covering the 
pseudo-rapidity ranges < 2.0 and |t]| < 1.4, respectively. The SPD contributes to the minimum-bias 
trigger if hits are detected on at least two different chips (each covering a detector area of 1 .28 x 1 .4 1 cm^) 
on the outer layer. The VZERO detector is composed of two anays of scintillator tiles covering the full 
azimuth in die pseudo-rapidity regions 2.8 < T] < 5.1 (VZERO-A) and -3.7 < T] < -1.7 (VZERO-C). 
The events used in this analysis were collected with two different interaction trigger configurations: in the 
first part of the data taking period, signals in two out of the three triggering detectors (SPD, VZERO-A, 
VZERO-C) were required, while in the second part a coincidence between the VZERO-A and VZERO-C 
detectors was used. Events were further selected offline to remove background coming from parasitic 
beam interactions on the basis of the timing information provided by the VZERO and the neutron ZDC 
detectors (two calorimeters located at z ±114 m from the interaction point). It was verified that the 
timing information from the ZDCs was available in all the hadronic interactions that passed the trigger 
condition. The luminous region had an r.m.s. width of about 6 cm in the longitudinal direction and 50- 
60 /im in the transverse direction. These values were stable during the entire data taking period. Only 
events with a vertex found within ±10 cm from the centre of the detector along the beam line were 
considered for the D meson signal extraction. 

Collisions were classified according to their centrality, defined in terms of percentiles of the hadronic 
Pb-Pb cross section and detemiined from the distribution of the summed amplitudes in the VZERO 
scintillator tiles. To obtain the total hadronic cross section, this distribution was fitted using the Glauber 
model for the geometrical description of the nuclear collision |[T5l complemented by a two-component 
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Table 1: Average values of the number of participating nucleons, and of the nuclear overlap function for the 
considered centrality classes, expressed as percentiles of the hadronic cross section. The values were obtained 
with a Monte Carlo implementation of the Glauber model assuming an inelastic nucleon-nucleon cross section of 
64 mb 121. 



Centrality class 


(A^part) 


{Taa) (mb-i) 


0-20% 


308 ±3 


18.93 ±0.74 


40-80% 


46±2 


1.20 ±0.07 


0-10% 


357 ±4 


23.48 ±0.97 


10-20% 


261 ±4 


14.43 ±0.57 


20-40% 


157 ±3 


6.85 ±0.28 


40-60% 


69±2 


2.00 ±0.11 


60-80% 


23 ±1 


0.42 ±0.03 



model for particle production 1129 11301 . The fit was performed in a range of measured VZERO amplitudes 
where the trigger is fully efficient for hadronic interactions and the contamination by electromagnetic 
processes is negligible Il30l . This range corresponds to 90±1% of the total hadronic cross section. The 
nuclear modification factor Raa was measured for D'', D+, and D*+ mesons as a function of transverse 
momentum for the centrality classes 0-20% and 40-80%. In order to study in more detail its central- 
ity dependence, /?aa was also evaluated, for wide pi intervals, in naixower centrality classes: 0-10%, 
10-20%, 20-40%, 40-60%, and 60-80%. Table [U shows the average values of the number of partic- 
ipating nucleons (A'pait), and of the nucleai" overlap function {Taa) in these centrality classes. In the 
centrality range considered in this analysis, 0-80%, and for both the configurations of the interaction 
trigger described above, the trigger and event selection are fully efficient for hadronic interactions, and 
the contamination by electromagnetic processes is negligible. 

In total, 13 X 10^ Pb-Pb collisions with centrality in the range 0-80% passed the selection criteria 
described above and were used in the analysis. The corresponding integrated luminosity is Lint = 
2.12±0.07 lJ.h-^. 

The trajectories of the D meson decay particles were reconstructed from their hits in the TPC and in 
the ITS. The TPC lISTTl provides track reconstruction with up to 159 three-dimensional space points per 
track in a cylindrical active volume that covers the region 85 < r < 247 cm and —250 < z < ±250 cm 
in the radial and longitudinal directions, respectively. The ITS |[32l consists of six cylindrical layers 
of silicon detectors with radii in the range between 3.9 cm and 43.0 cm. Around the two innermost 
layers equipped with pixel detectors (SPD, described above), Silicon Drift Detectors (SDD) are used 
in the two intermediate layers, while the two outermost layers are made of double-sided Silicon Strip 
Detectors (SSD). The alignment of the ITS sensor modules, which is crucial to achieve the high space 
point resolution needed in heavy flavour analysis, was performed using survey information, cosmic-ray 
tracks, and pp data, with the methods described in ||32| . For the residual misalignment along the 
coordinate, an r.m.s. of about 8 /im for SPD and 15 /im for SSD modules was estimated Il32ll33l . For 
SDD, with the calibration level reached on the 2010 data sample, the space point resolution along is 
60 nm for those modules that do not suffer from significant drift field non-uniformities. The residual 
misalignment is included in an effective way in the detector simulation by randomly displacing the ITS 
modules with respect to their ideal positions according to the estimated precision of the alignment. 

The primary vertex position and covariance matrix were determined from the tracks reconstructed in the 
TPC and ITS by using an analytic minimization method, applied after approximating the tracks to 
straight lines in the vicinity of their common origin. The same algorithm was used for the reconstruction 
of the decay vertices of D'' and D+ candidates. The high spatial resolution of the reconstmcted hits, 
together with the low material budget (on average 7.7% of a radiation length for the ITS at T] = 0) and 
the small distance of the innermost layer from the beam vacuum tube, allows for the measurement of the 
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Figure 1: Left: track impact parameter resolution in the transverse plane as a function of pt in Pb-Pb collisions. 
For Pt <2 GeV /c, pion identification by the TPC or TOF detectors is required; the results for data and simulation 
are shown. The simulation includes the effect of the residual geometrical misalignment of the sensor modules of 
the Inner Tracking System (see text for details). Right (colour online): difference between the measured time- 
of-flight and that expected under the kaon hypothesis as a function of track momentum for the 20% most central 
Pb-Pb collisions. 



track impact parameter in the transverse plane (do), i.e. the distance of closest approach of the track to 
the primary vertex along r(j), with a resolution better than 65 /im for transverse momenta pt> I GeV/c. 
The impact parameter resolution Od^ is shown in the left-hand panel of Fig.[T]as a function of pi, for data 
and simulation, for charged hadron tracks selected with the same criteria used in the D meson analysis. 
The applied track quality cuts were based on the request of having at least 70 associated space points 
(out of a maximum of 159) in the TPC with a per degree-of-freedom of the momentum fit lower than 
2, and at least 2 associated hits in the ITS, out of which at least one has to be in the silicon pixel layers. 
Only tracks with ti'ansverse momentum pt > 0.5 GeV/c (0.7 for the 20% most central collisions) and 
|t] I < 0.8 were used for the D meson analysis and are displayed in Fig. [T](left). For pt <2 GeV/c, only 
particles identified as pions were selected, as explained in the next paragraph. The impact parameter 
resolution is better than for pp collisions lETl . e.g. by 10 /im at pt = I GeV/c, since, in the Pb-Pb 
case, the primary vertex is reconstructed using a larger number of tracks, hence with better precision. 
Indeed, the resolution on the transverse coordinates of the primary vertex is about 5 /im in central Pb-Pb 
collisions, while it is up to 40 /im in pp ll27l . The systematic effect on the D meson analysis of the 
small difference in resolution (5 /im) between data and simulation will be discussed in Section [6l The 
resolution on the transverse momentum of tracks reconstructed in the TPC and in the ITS, and passing 
the quality selection criteria described above, was measured to be about 1% at ;?t = 1 GeV/c and about 
2% at pt = lO GeV/c. 

The particle identification (PID) capabilities are provided by the measurement of the specific energy 
loss dE/dx in the TPC and of the time-of-flight in the TOF detector. The dE/dx samples measured 
by the TPC aie reduced, by means of a truncated mean, to a Gaussian distribution with a resolution of 
(^dE lAxI (d^" / dx) ~ 6% which is shghtly dependent on track quality and detector occupancy. The TOF 
detector 041 is positioned at 370-399 cm from the beam axis and covers the full azimuth and the pseudo- 
rapidity range \v\\< 0.9. In Pb-Pb collisions, in the centrality range 0-70%, the overall time-of-flight 
resolution was measured to be about 90 ps for pions with a momentum of 1 GeV/c. This value includes 
the detector intrinsic resolution, the electronics and calibration contribution, the uncertainty on the start 
time of the event, and the tracking and momentum resolution. The start time of the event is measured by 
the TO detector, made of two arrays of Cherenkov counters located on either side of the interaction point 
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and covering the pseudorapidity ranges —3.28 < rj < —2.97 and 4.61 < rj < 4.92, respectively. For the 
events in which the TO signal is not present, the start time is estimated using the particle arrival times 
at the TOF. In the centrality class 70-80%, the TOF resolution slightly worsens due to the increasing 
uncertainty on the start time determination, while still remaining below 100 ps. 

In the right-hand panel of Fig. [T] the difference between the measured time-of-flight and that expected 
under the kaon hypothesis is shown as a function of the track momentum for Pb-Pb collisions in the 
centrality range 0-20%, illustrating the separation between pions, kaons and protons. The expected 
time-of-flight is calculated for the given mass hypothesis from the total integrated path length and the 
measured momentum of the track. A compatibility cut with the PID response from the TPC was used 
in order to decrease the contamination from tracks with wrong hit association in the TOF detector. The 
bands corresponding to particles identified as pions (upper band), kaons (middle) and protons (lower 
band) are separated up to p « 2 GeV/c, corresponding to the momentum range in which the TOF PID is 
used in this analysis. 



3 D meson reconstruction and selection 

The D°, D+, and D*+ mesons and their antiparticles were reconstructed in the central rapidity region 
from their charged hadronic decay channels D° — K^;r+ (with branching ratio, BR, of 3.87 ±0.05% 
and mean proper decay length cT « 123 /im), D+ — )• K"7r+7r+ (BR of 9.13 ±0.19%, cT w 312 /im), 
and D*+ — )• D°7r+ (strong decay with BR of 67.7 ±0.5%) Il35l . The D meson yields were extracted from 
an invariant mass analysis of fully reconstructed decay topologies displaced with respect to the primary 
vertex, using the same procedure as for pp collisions lETI . 

D° and D+ candidates were defined from pairs and triplets of tracks with proper charge sign combination 
and selected by requiring at least 70 associated space points in the TPC, with ;t;^/ndf < 2, and at least 
2 associated hits in the ITS, out of which at least one in the SPD. A fiducial acceptance cut < 0.8 
was applied as well, along with a transverse momentum threshold pi > 0.5 GeV/c (0.7 for the 20% most 
central collisions), aimed at reducing the large combinatorial background. 

D*+ candidates were obtained by combining the D" candidates with tracks selected with transverse mo- 
mentum pi > 0.2 GeV/c in the centrality range 0-20% and pt>0.l GeV/c in 20-80%. The momentum 
of the pion from the D*+ decay is typically low, because of the small mass difference between the D*+ 
and D" mesons. In order to reduce the combinatorics, these tracks were selected requiring at least 3 
associated ITS hits (4 in the 0-20% centrality class), in addition to the same TPC quality selection as 
that used for the D" and decay tracks. In the centrality class 40-80%, also tracks reconstructed only 
in the ITS, with at least 3 hits, were used to enhance the D*+ signal at low pi. 

The selection of the D° and D+ decays was based on the reconstruction of secondary vertex topologies, 
with a separation of a few hundred microns from the interaction point. In the case of the D*+ decay, 
the secondary vertex topology of the produced D° was reconstructed. The selection is essentially the 
same as that used for the pp case ETl and exploits the separation between the secondary and primary 
vertices (decay length) and the pointing of the reconstructed meson momentum to the primary vertex. 
The pointing condition is applied by requiring a small value for the angle ^pointing between the directions 
of the reconstructed momentum of the candidate and of its flight line, defined by the positions of the 
primary and secondary vertices. In order to cope with the much larger combinatorial background and 
to exploit the better resolution on the reconstructed primary vertex position, the cuts were in general 
tightened with respect to the pp case. Two additional cuts, on the projections of the pointing angle and 
of the decay length in the transverse plane (0poj„ting ^^'^ ^^)' were introduced to further suppress the 
combinatorial background. 

The cuts were defined so as to have large statistical significance of the signal and to keep the selection 
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efficiency as high as possible. This latter requirement was dictated also by the fact that too tight cuts 
result in an increased contribution to the raw yield from feed-down D mesons originating from decays 
of B mesons. It was also checked that background fluctuations were not causing a distortion in the 
signal line shape by verifying that the D meson mass and its resolution were in agreement with the PDG 
value and the simulation results, respectively. The resulting cut values depend on the D meson pi and 
on the centrality of the event. They lead to a selection efficiency that increases with increasing pi and 
decreases from peripheral to central collisions: looser cuts could be used for peripheral events, where the 
combinatorial background is lower. The cut values quoted in the following refer to the tightest selections 
in the lower pi intei^vals for the 0-20% centrality class. 

The PID selection relies on the pion and kaon identification by the TPC and TOP detectors. Cuts at ±3 a 
around the expected mean energy deposit AE /dx and time-of-flight were used. This selection provides 
a strong reduction, by a factor of about 3, of the combinatorial background in the low-pt region, while 
preserving most of the signal (ss 95% according to simulations, as detailed in the next Section). In 
the D*+ case, a tighter PID cut at 2 a on the TPC dE/Ax was applied to the D" decay products in the 
centrality class 0-20%, in order to cope with the large combinatorial background. 

With the track selection described above, the acceptance in rapidity for D mesons drops steeply to zero 
for \y\ > 0.5 at low pi and \y\ > 0.8 for px^5 GeV/c. A p^ dependent fiducial acceptance cut was therefore 
applied on the D meson rapidity, \y\ < ymiPt), with 3'fid(pt) increasing from 0.5 to 0.8 in < ;?t < 
5 GeV/c according to a second order polynomial function and taking a constant value of 0.8 for pi > 
5 GeV/c. 

For D" mesons, the two decay tracks were selected requiring a significance of the impact parameter 
with respect to the event primary vertex |J()|/arf„ > 0.5 and a maximum distance of closest approach 
between each other of 250 /^m. The minimum decay length was set at 100 /im. Furthermore, the cuts 
X < —45000 /im^ on the product of the decay track impact parameters and L^/an? > 7 on the 
significance of the projection of the decay length in the transverse plane were applied. A selection on the 
angle 6* between the kaon momentum in the D*' rest frame and the boost direction was used to reduce 
the contamination of background candidates that do not represent real two-body decays and typically 
have large values of | cos 6* \ . The applied cut was | cos d* \ < 0.8. The pointing of the D° momentum to 
the primary vertex was imposed via the cuts cos ^pointing > 0.95 and cos ^p^jntj^g > 0.998. 

For D+ mesons, a decay length of at least 1.9 mm was required. It should be noted that D+ mesons in 
the 0-20% centrality class are reconstructed only for pt > 6 GeV/c, where the Lorentz dilation of the 
D+ lifetime allows for a tight cut on the decay length. Further requirements to reduce the combinatorial 
background were cos 0pointing > 0.995, cosd-^^-^^^^ > 0.997, V^' /Ot^y > 12, and Y^d^ > (300 /im)^ (sum 
of the squared impact parameters of the three decay tracks). The D+ cuts ai^e in general tighter than the 
D" ones because of the larger combinatorial background. 

In the D*+ analysis, the selection of the decay D° was similar to that used for the D" analysis, with a 
tighter cut on the pointing angle, cos ^pointing > 0.99. The decay pion was selected with the track quality 
cuts described above and requiring a minimum pt that varied in the range 0.1-1 GeV/c depending on the 
D*+ momentum and event centrality. In the 0-20% centrality class and for D*+ transverse momentum 
below 6 GeV/c, a 3a compatibility cut with respect to the pion expectation values was applied to the 
measured dE/dx and time-of-flight. 

Figure |2] shows the invariant mass distributions of the selected D*', D+, and D*+ candidates in some of 
the pi intervals used in the analysis, for the 0-20% centrality class. The D*' and D+ yields were extracted 
by fitting the distributions with a function composed of a Gaussian for the signal and an exponential 
term that describes the background shape. The D*+ background was described with a threshold function 
multiplied by an exponential ETl . The centroids of the Gaussians were found to be compatible with 
the PDG masses of the D mesons (35l, and their widths to be well reproduced in the simulation. The 
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Figure 2: Invariant mass distributions for (upper panels), D+ (central panels), and D*+ (lower panels) candi- 
dates and their charge conjugates in selected pt intervals for 3.2 x 10^ 0-20% central Pb-Pb collisions. The curves 
show the fit functions described in the text. The values of mean (ji) and width (a) of the signal peak are reported 
in the plots together with the raw signal yield. The uncertainties on the signal yields reported in the figures are 
statistical only. 



signal yields (sum of particle and antiparticle) are reported in Table |2] for the pt intervals considered in 
the analysis, for the centrality classes 0-20% and 40-80%. 



4 Corrections 

The D meson raw yields extracted from the fits to the invariant mass distributions were corrected to 
obtain the production yields for primary (i.e. not coming from weak decays of B mesons) D°, D+, and 
D*+. The contribution of secondary D mesons from B decays, which is of the order of 15% as explained 
in the following, was estimated using pQCD predictions for B production and Monte Carlo simulations. 
The D mesons remaining after the subtraction of the B feed-down contribution are those produced at the 
interaction vertex, and they will be referred to as 'prompt' in the following. 
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Table 2: Measured raw yields for D", D+, and D*+ mesons and their antiparticles in the transverse momen- 
tum intervals considered for the 0-20% and 40-80% centrality classes. The systematic uncertainty estimation is 
described in Section |6] 



interval 
(GeV/c) 





-20% centrality 
D++D 


A?''^* ±St£ 

D*++D* 


It. ± syst. 

4 

D0 + D° 


0-80% centrality 
D+ + D 


D*++D* 


2-3 


538 ± 84 ±43 






231±31±12 




82±21±12 


3-4 


774 ±108 ±46 






241 ±32 ±12 


58±19± 9 


4-5 


583 ± 79 ±35 




60±18±12 


176±20± 9 


114±22± 6 


36± 7± 5 


5-6 


318± 67±19 




87±13± 4 


29± 9± 3 


6-8 


342 ± 48 ±21 


167 ±43 ±33 


63±16± 6 


113±14± 6 


130 ±34 ±20 


47±13± 5 


8-12 


327 ± 41 ±20 


132 ±38 ±20 


55±12± 6 


107±15± 6 


119±26±18 


57±11± 6 


12-16 


67 ± 15 ± 7 


62±15± 6 


38 ± 8± 4 


41± 9± 2 




23 ± 6± 2 



The prompt D meson production yields were calculated starting from the raw yields (N^^"", reported 
in the previous section) divided by a factor of two to evaluate the charge (paiticle and antipaiticle) 
averaged yields. These were corrected for the B meson decay feed-down contribution (i.e. multiplied 
by the prompt fraction /prompt)> and divided by the acceptance-times-efliciency for prompt D mesons, 
(Acc X £)prompt- They were normalized according to the decay channel branching ratio (BR), pt interval 
width (Apt), rapidity coverage (Ay), and the number of events analyzed (A^evt)- As an illustration, the D+ 
yields were computed as: 

1 /prompt(A)-5A^°*""(A) , 
AyApt (AcCX£)prompt(A)-BR-A^evt 

As mentioned in Section[3j the D meson yields were measured in a rapidity range varying from \y\ < 0.5 
at low Pt to \y\ < 0.8 at high pt. The rapidity acceptance correction factor Ay = 2yM assumes a uniform 
rapidity distribution for D mesons in the measured y range. This assumption was checked to the 1% 
level lizl with PYTHIA IMl pp simulations with the Pemgia-0 tuning [37 1. 

The acceptance-times-efficiency corrections Acc x e were obtained using Monte Carlo simulations. 
Minimum-bias Pb-Pb collisions at ^/s^ = 2.76 TeV were produced with the HIJING vl.36 event gen- 
erator Il38l . Prompt and feed-down (B decays) D meson signals were added using pp events from the 
PYTHIA v6.4.21 event generator [36] with the Perugia-0 tuning [1371 . Each injected pp event was re- 
quired to contain a cc or bb pair and D mesons were forced to decay in the hadronic channels of interest 
for the analysis. The number of pp events added to each Pb-Pb event was adjusted according to the 
Pb-Pb collision centrality. The simulations used the GEANT3 |[39l particle transport package together 
with a detailed description of the geometry of the apparatus and of the detector response. The simulation 
was configured to reproduce the conditions of the luminous region and of all the ALICE subsystems, 
in terms of active electronic channels, calibration level, and their time evolution within the Pb-Pb data 
taking period. 

The efficiencies were evaluated in centrality classes corresponding to those used in the analysis of the data 
in terms of charged-particle multiplicity, hence of detector occupancy. Figure |3] shows the D" — )• K^;r+, 
D"*" — K^7r+7r+, and D*+ — D'';r+ acceptance-times-efficiency for prompt and feed-down D mesons 
with rapidity \y\ < yM- The efficiencies correspond to Pb-Pb collisions in the centrality class 0-20%. 
The selection cuts described in Section |3] were applied. The values for the case of not applying PID 
aie shown as well, in order to point out that this selection is about 95% efficient for the signal. For the 
three meson species, the acceptance-times-efficiency increases with pt, starting from few per mil and 
reaching 5-10% at high pt. No significant difference in the acceptance-times-efficiency for particles 
and antiparticles was observed. 



dpt 



Suppression of high pt D mesons in Pb-Pb collisions at ^^s^ = 2.76 TeV 
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Figure 3: Acceptance-times-efficiency in Pb-Pb collisions (0-20% centrality class) for (left), D+ (middle), and 
D*+ (right) mesons. The efficiencies for prompt (sohd lines) and feed-down (dotted fines) D mesons are shown. 
Also displayed, for comparison, the efficiency for prompt D mesons without PID selection (dashed lines). 

The acceptance-times-efficiencies for D mesons from B decays are larger than for prompt D mesons by 
a factor of approximately 2, because the decay vertices of the feed-down D mesons are more displaced 
from the primary vertex and, thus, more efficiently selected by the cuts. 

In the 40-80% centrality class, as discussed in the previous Section, the selection cuts were looser, 
resulting in a higher efficiency. The dependence of the D meson selection efficiency on the detector oc- 
cupancy was evaluated by comparing the efficiencies for central (0-20% centrality class) and peripheral 
(40-80% centrality class) events when applying the same selection cuts (those of the 0-20% class were 
used as a reference). The results showed only small variations as a function of centrality, e.g. 5-10% 
for D°, as expected from the small variation of the single track reconstruction efficiency with central- 
ity ifTSl . Indeed, also the efficiency of the topological selection is expected to be practically independent 
of centrality in the considered range 0-80% where the resolution on the primary vertex position is not 
significantly affected by the multiplicity of tracks used in its determination. 

The prompt D meson production yields dN/dpt in Pb-Pb collisions were obtained by subtracting the 
contribution of D mesons from B decays with the same procedure used for the measurement of the 
production cross sections in pp collisions lITTl . In detail, the feed-down contribution was estimated 
using the beauty production cross section from the FONLL calculation BOl . the B— >D decay kinematics 
from the EvtGen package PTI . and the Monte Carlo efficiencies for feed-down D mesons. For Pb- 
Pb collisions, the FONLL feed-down cross section in pp at ^/s = 2.76 TeV was scaled by the average 
nuclear overlap function (Taa) in each centrality class. Thus, omitting for brevity the symbol of the 
7?t-dependence (pt), the fraction of prompt D mesons reads: 

/p,„^pt = 1-{N^ feed-down raw/^D raw) ^ 

/ d^q \ „feed-down (AcC X £)feed-down " Apt • BR • A^evt 

VdjdAA-eed-down" ' ^^^2 



(Taa) 



where (Acc x £)feed-down is the acceptance-times-efficiency for feed-down D mesons. The nuclear mod- 
ification factor of the feed-down D mesons, 7?^^''^'^°*", is related to the nuclear modification of beauty 
production in Pb-Pb collisions, which is cuixently unknown. We therefore assumed for the coixection that 
the nuclear modification factors for feed-down and prompt D mesons are equal (7?^^''^'^°^" = 7?^™'") and 
varied this hypothesis in the range 1/3 < 7?^^'*^'*°*'"//?^"'"^' < 3 to deteiTnine the systematic uncertainty. 
This hypothesis is justified by the range of the model predictions for the charm and beauty Raa II101I16II 
and, as discussed in Section [6l by the CMS Collaboration results on 7?aa for non-prompt J/y |!25|. The 
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value of /prompt depends on the D meson species, the transverse momentum intei^val, the applied cuts, 
the parameters used in the FONLL B prediction, and the hypothesis on 7?^^''^''°*". The resulting val- 
ues, for the case 7?^^''^'^°^" = 7?^™'", range from 0.95 in the lowest transverse momentum interval 
(2 < ;7t < 3 GeV/c) to « 0.85 at high pf 



5 Reference pp cross section at ^ = 2.76 TeV 

The reference pp cross sections used for the determination of the nuclear modification factors were ob- 
tained by applying a -y/i-scaling H2l to the cross sections measured at ^/s = 7 TeV ETIl . The scaling 
factor for each D meson species was defined as the ratio of the cross sections from the FONLL pQCD 
calculations |@Q at 2.76 and 7 TeV. The same values of the pQCD factorization scale /Xp and renormal- 
ization scale /Xr, and of the charm quark mass nic were used in the calculation for the different energies. 
Namely, jLtp = Mr = "it with nii = y'pf+rn^ and Mc = 1.5 GeV/c^. The theoretical uncertainty on the 
scaling factor was evaluated by considering the envelope of the scaling factors resulting by varying in- 
dependently the scales in the ranges 0.5 < /iR/mt < 2, 0.5 < /^p/'Mt < 2, with 0.5 < Hr/Hf < 2, and the 
quark mass in the range 1.3 < < 1.7 GeV/c^, following the prescription in B3l . This uncertainty 
ranges from 1 jo% at pt = 2 GeV/c to about ±5% for pt > 10 GeV/c |42l . The procedure was validated 
by scahng the ALICE pp data to the Tevatron energy, ^/s = 1.96 TeV, and comparing to CDF measure- 
ments II42[|44II . In addition, it was verified that the scaling factor and its uncertainty are the same if the 
GM-VFNS calculation El is used instead of FONLL Il42l . 

The D", D+, and D*+ cross sections were measured, though with limited precision and pt coverage, in 
pp collisions at ^/s = 2.76 TeV using a sample of about 6 x 10^ minimum-bias events collected during 
a short run at the same energy as Pb-Pb collisions. These measurements were found to be in agreement 
with the scaled 7 TeV measurements, within statistical uncertainties of about 20-40% depending on pi 
and on the meson species |[28l . 



6 Systematic uncertainties 

Systematic uncertainties on the Pb-Pb yields 

The systematic uncertainties on the prompt D meson yields in Pb-Pb collisions are summarized in Table |3] 
for the lowest and highest p^ intervals in the two centrality classes 0-20% and 40-80%. 

The systematic uncertainty on the yield extraction from the invariant mass spectra was determined by 
repeating the fit, in each p^ interval, in a different mass range and also with a different function to 
describe the background. Namely, a parabola, instead of an exponential, was considered for D*' and D+, 
and a power law multiplied by an exponential or a polynomial for D*+. A method based on counting the 
signal in the invariant mass distribution, after subtraction of the background estimated from a fit to the 
side bands, was also used. The uncertainty was defined as the maximum difference of these results and 
it was found to vary in the range 5-20%, depending on the pt interval and on the collision centrality. 

The systematic uncertainty on the tracking efficiency was estimated by comparing the efficiency (i) of 
track finding in the TPC and (ii) of track prolongation from the TPC to the ITS between data and simula- 
tion, and (iii) by vaiying the track quality selections. The efficiency of track prolongation from the TPC 
to the ITS and of association of hits in the silicon pixel layers was found to be described in simulation 
at the level of 5% in the pt range relevant for this analysis (0.5-15 GeV/c). The centrahty dependence 
of these efficiencies, which is limited to ±3% in this pt range, was found to be reproduced within 1.5%. 
The effect of wrong association of ITS hits to tracks was studied in the simulation. It was found that the 
fraction of D mesons with at least one decay track with wrong hit associations increases with centrality, 
due to the higher detector occupancy, and vanishes at large pt, where the track extrapolation between lay- 
ers is more precise. In the centrality class 0-20%, it ranges from 7% to 1% in the transverse momentum 
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Table 3: Summary of relative systematic uncertainties on the prompt D meson production yields in Pb-Pb colli- 
sions for the lowest and highest pt bins measured for the three mesons. 



Particle 







D+ 


D*+ 






z, — J) 




6-8 12-16 


4-6 


12-16 




Yield extraction 


8% 


10% 


20% 


10% 


20% 


10% 




Tracking efficiency 


10% 


10% 


15% 


15% 


15% 


15% 


0-20% 
centrality 


rVD cincicncy 


13% 

^>^3% 

4% 


10% 

5% 
3% 


15% 

5% 
1% 


15% 

5% 
5% 


10% 

3% 


10% 

5% 
3% 




FONLL feed-down com 


+ 2% 

-14 


-1% 


+?% 




ti% 






y^feed-down /y^prompt 


+ '^% 
-10 




+ ■^% 

-16 




tt% 






RT? 
D K 


1.3% 


2.1% 


1.5% 




v_-CllLld.ilLy iillllLS 


<0.1% 




iTitf*r\/'al ( (~\p\l 1 r'\ 

IJ^ lllLti \ (XL v,vJt y 1 L. J 


2-3 


12-16 


3-4 


8-12 


2-4 


12-16 




iieiQ exiracuon 


5% 


5% 


15% 


15% 


15% 


8% 




irdCKjng eiiiciency 


10% 


10% 


15% 


15% 


15% 


15% 


40-80% 
centrality 


Cut efficiency 
PID efficiency 


13% 


10% 

5% 


10% 

5% 


10% 

5% 


10% 

;io% 


10% 
5% 




MC Pt shape 


1% 


3% 


1% 


3% 


5% 


4% 




FONLL feed-down corr. 


+ 3% 

-16 


^> 




tt% 


1 1% 


11% 




^feed-down/^prompt 


+ 5% 
-12 






+ ''% 

-20 


tl% 


1^% 




BR 


1.3% 


2.1% 


1.5% 




Centrality limits 


3% 



interval 2<pi<\6 GeV/c. However, it was verified that the signal selection efficiencies are compatible, 
within statistical uncertainties, between D mesons with and without wrong hit associations. Indeed, the 
mis-associated hit is typically very close in space to the coixect hit. Overall, the systematic uncertainty 
from track reconstruction amounts to 5% for single tracks, which results in a 10% uncertainty for D*' 
mesons (two-track final state) and 15% uncertainty for D+ and D*+ mesons (three-track final state). 

The uncertainty on the correction for the selection cut efficiency was evaluated by repeating the analysis 
with different sets of cuts and was defined as the variation of the resulting corrected yields about the 
value corresponding to the central set. This resulted in 13% for D° for p\_ <3 GeV/c, 15% for D+ in 
all Pt intervals in the 0-20% centrality class, and 10% for the other cases (see Table |3]l. Part of this 
uncertainty comes from residual detector misalignment effects not fully described in the simulation. In 
order to estimate this contribution, the secondary vertices in the simulation were also reconstructed after 
scaling, for each track, the impact parameter residuals with respect to their true value. In particular, a 
scahng factor of 1.1-1.2 was applied in order to reproduce the impact parameter resolution observed in 
the data (see Fig.[T|). The relative variation of the efficiency is 8% for pt = 2-3 GeV/c and negligible for 
Pt > 5 GeV/c. This effect was not included explicitly in the systematic uncertainty, since it is already 
accounted for in the cut variation study. A further check was performed by comparing the distributions 
of the cut variables used for the candidate selection in the data and in the simulation. These comparisons 
can only be carried out by releasing the selection, hence essentially for background candidates. However, 
they provide an indication of the level of accuracy of the simulation. A good agreement was observed. 
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with no dependence on collision centrality. 

The uncertainty arising from the PID selection was estimated by comparing the corrected signals ex- 
tracted with and without this selection. In the 20% most central collisions, it was found to be t^l'^^ for 
Pt<6 GeV/c and ±5% for pt> 6 GeV/c. In the 40-80% centrality range, it was estimated to be ±5% 
for pt>3 GeV/c and t^^l m2<pt<3 GeV/c. 

The uncertainty on the efficiencies arising from the difference between the real and simulated shape of the 
D meson transverse momentum distribution, which includes also the effect of the pt dependence of the 
nuclear modification factor, depends on the width of the pi intei-vals and on the slope of the efficiencies 
with Pt. It was estimated by varying the simulated shape between the PYTHIA and FONLL dN/dpt, 
with and without the nuclear modification observed in the data. The resulting uncertainty is below 5% in 
all the Pt intervals considered for the three meson species. As an example, for D'' it is 4% in the lowest 
and highest pt intervals (2-3 GeV/c and 12-16 GeV/c) and 1% in 3-12 GeV/c. 

The 7?t-differential yields for D'' and mesons, extracted separately, were found to be in agreement 
within the statistical uncertainties of about 20-25%. Due to the limited statistics, this check could not be 
carried out for D+ and D*+ mesons. 

The systematic uncertainty from the subtraction of feed-down D mesons from B decays was estimated as 
for the pp case ll27l . The contribution of the FONLL perturbative uncertainties was included by varying 
the heavy quark masses and the factorization and renormaUzation scales in the ranges proposed in B3l . 
Furthermore, a different procedure was used to evaluate the prompt fraction. In this approach, the ratio 
of the FONLL feed-down and prompt production cross sections is the input for evaluating the coiTcction 
factor. Then, the prompt fraction depends explicitly on the ratio of nuclear modification factors of feed- 
down and prompt D mesons: 



prompt 



1 + 



V 



(ACC X £)feed-down 
(Acc X £)prompt f cPa 



N FONLL 

V^yd]^ J feed-down 
FONLL 

3^ 7 prompt 



R 



feed— down 
AA 



-1 



prompt 
AA 



(4) 



The systematic uncertainty due to the B feed-down subtraction was evaluated as the envelope of the 
results obtained with the two methods, Eqs. ^ and (|4]i, when vaiying the FONLL parameters. The 
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Figure 4: Relative variation of the prompt D'' meson yield as a function of the hypothesis on R^^^ down^^prompt 
for the B feed-down subtraction approach based on Eq. (|3). 
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Table 4: Summary of relative systematic uncertainties 
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on Raa- For the data systematic uncertainties and the B 
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D+ 
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resulting uncertainty ranges between 't[4% at low pt and lg% at high pt (for D° in the 0-20% centrality 
class). 

The contribution from the different nuclear modification factors of prompt and feed-down D mesons was 
evaluated by varying the hypothesis on 7?^^''-'^°^7/?p"'"p' in the range 1 /3 < < 3 for 

both feed-down subtraction methods. The resulting uncertainty is at most 30%, as shown in Fig.|4l where 
the relative prompt D" yield variation is displayed as a function of /?^5f-'*°*V/?P;°"'P' for four pt intervals 
using the B feed-down subtraction approach based on Eq. Q. Considering the resulting values of 7?^™'" 



shown in the next section, the variation of the hypothesis on ^^a^ 



/^AA™'" corresponds, for the 20% 



most central collisions, to values of the nuclear modification factor of D mesons from B feed-down in a 
range of about 0.17-1.5 at low pi and 0.09-0.8 at high pt. The /?aa of non-prompt J/y^ measured by 
CMS |[25l . falls in this range as well as the available model predictions for B meson energy loss II10II161 . 

The contribution due to the 1.1% relative uncertainty on the fraction of hadronic cross section used in the 
Glauber fit to determine the centrality classes was obtained by estimating the variation of the D meson 
dN/dpt when the limits of the centrality classes are shifted by ±1.1% (e.g., for the 40-80% class, 40.4— 
80.9% and 39.6-79.1%). The resulting uncertainty is common to all meson species and all pi bins for a 
given centrality class. It increases from central to peripheral events. In particular, it is less than 0.1% in 
the 0-20% centrality class and 3% in 40-80%. 

Finally, the systematic uncertainty on the branching ratios |[35l was considered. 
Systematic uncertainties on /?aa 



The systematic uncertainties on the /?aa measurement derive from the uncertainties on: the reference 
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cross section for pp collisions, the Pb-Pb yields, and the average nuclear overlap function for the various 
centrality classes, as given in Table [T] For the pp reference, the uncertainties on the measurement at 
= 7 TeV were quantified in ETl and the scaling to = 2.76 TeV, described in Section[5J introduces 
additional uncertainties of about 10-30%. The uncertainties on the Pb-Pb prompt D meson yields were 
described previously in this section. For the nuclear modification factor, the pp and Pb-Pb uncertainties 
were added in quadrature, except for the feed-down contribution deriving from FONLL uncertainties, 
that partly cancels in the ratio. This contribution was evaluated by comparing the Raa values obtained 
with the two methods for feed-down correction described above and with the different heavy quark 
masses, factorization and renormalization scales used in FONLL. In this study, the same method and the 
same set of FONLL parameters were used for pp and Pb-Pb, so as to take into account the correlations 
of these sources in the numerator and denominator of Raa- 

The resulting systematic uncertainties are summarized in Table |4] In the table, the normalization uncer- 
tainty is the quadratic sum of the 3.5% pp normalization uncertainty ETl . the contribution due to the 
1.1% uncertainty on the fraction of hadronic cross section used in the Glauber fit discussed above, and 
the uncertainty on (Taa), which is 3.9% for the centrality class 0-20% and 5.9% for the 40-80% class. 

7 Results 

7.1 D meson pt spectra and Raa 

The transverse momentum distributions dN/dpt of prompt D^, D+, and D*+ mesons are presented in 
Fig. m for the centrality classes 0-20% and 40-80%. The spectra from Pb-Pb collisions, defined as 
the feed-down corrected production yields per event (see Eq. Q), are compared to the reference spectra 
from pp collisions, which are constructed as (Taa) da/dpt, using the y(s-scaled pp measurements at 
7 TeV ETl and the average nuclear- overlap function values from Table[T] A clear suppression is observed 
in Pb-Pb collisions, which is stronger in central than in peripheral collisions. 

The ratio of the Pb-Pb to the reference spectra provides the nuclear modification factors RAA{Pt) of 
prompt D°, D+, and D*+ mesons, which are shown for central (0-20%) and semi-peiipheral (40-80%) 
collisions in Fig. |6l The vertical bars represent the statistical uncertainties, typically about 20-25% for 
D° and about 30-40% for D+ and D*+ mesons in central collisions. The total pt-dependent systematic 
uncertainties, shown as empty boxes, include all the contributions described in the previous section, 
except for the normalization uncertainty, which is displayed as a filled box at Raa = 1 ■ The results for 
the three D meson species are in agreement within statistical uncertainties and they show a suppression 
reaching a factor 3^ {Raa ~ 0.25-0.3) in central collisions for pt > 5 GeV/c. For decreasing pi, the D'' 
Raa in central collisions shows a tendency to less suppression. 

The centrality dependence of the nuclear modification factor was studied in the two wider transverse 
momentum intervals 2 < pt <5 GeV/c, for D'^, and 6 < pt < 12 GeV/c, for the three D meson species. 
This study was performed in five centrality classes from 0-10% to 60-80% (see Table [T|). The invariant 
mass analysis and all the corrections were carried out as described in Sections |3] and H] The systematic 
uncertainties are essentially the same as for the pt-dependence analysis, except for the contribution from 
the D meson pt-shape in the simulation, which is larger in the wide intervals. It amounts to 8% for 
D°, 10% for D+, and 5-15% (depending on centrahty) for D*+ mesons in 6 < pt < 12 GeV/c. In the 
transverse momentum interval 2-5 GeV/c, this uncertainty is larger (8-17%, depending on centrality) 
due to the larger contribution from the pi dependence of the nuclear modification factor. The resulting 
Raa is shown in Fig.|7]as a function of the average number of participants, {Npart}- The contribution to the 
systematic uncertainty that is fully correlated between centrality classes (normalization and pp reference 
cross-section) and the remaining, uncorrected, systematic uncertainties are displayed separately, by the 
filled and empty boxes, respectively. The contribution from feed-down correction was considered among 
the uncorrected sources because it is dominated by the variation of the ratio 7?^^'*^''°*"//?^'"^^', which 
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Figure 5: (colour online) Transverse momentum distributions dN/dpt of prompt D" (left) and D+ (centre), and 
D*+ (right) mesons in the 0—20% and 40—80% centrality classes in Pb-Pb collisions at -y^iNN = 2.76 TeV. The 
reference pp distributions (Taa) da/dpt are shown as well. Statistical uncertainties (bars) and systematic uncer- 
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may depend on centrality. For the px interval 6-12 GeV/c, the suppression increases with increasing 
centrality. It is interesting to note that the suppression of prompt D mesons at central rapidity and high 
transverse momentum, shown in the right-hand panel of Fig. |7]is very similar, both in size and centiality 
dependence, to that of prompt J/\j/ mesons in a similar pt range and \y\ < 2.4, recently measured by the 
CMS Collaboration flBl. 
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7.2 Comparisons to light-flavour hadrons and with models 

In this section, the average nuclear modification factor of the three D meson species is compared to that 
of charged particles |[T9l . mainly light-flavour hadrons, and to model calculations. The contributions of 
D°, D+, and D*+ to the average were weighted by their statistical uncertainties. Therefore, the resulting 
Raa is close to that of the D" meson, which has the smallest uncertainties. The systematic errors were 
calculated by propagating the uncertainties through the weighted average, where the contributions from 
the tracking efficiency, from the B feed-down connection, and from the FONLL scaling of 7 TeV data 
to 2.76 TeV were taken as fully correlated among the three D meson species. The possible statistical 
correlation between the D" and D*+ Raa, induced by the D*+ — ^ D°;r+ decay, is negligible, because the 
statistical uncertainties, used as weights, are mainly determined by the background uncertainties, which 
are uncorrelated. The resulting values are shown in Table |5] for the two centrality classes where Raa 
was measured as a function of pt, and in Table [6] for the Raa as a function of centrality in the transverse 
momentum range 6 < pt < 12 GeV/c. 

In addition to final state effects, where parton energy loss would be predominant, also initial-state effects 
are expected to influence the measured Raa- In particular; the nuclear modification of the parton distribu- 
tion functions of the nucleons in the two colliding nuclei modifies the initial hard scattering probability 
and, thus, the production yields of hard partons, including heavy quarks. In the kinematic range relevant 
for chaim production at LHC energies, the main expected effect is nuclear- shadowing, which reduces 
the parton distribution functions for partons with nucleon momentum fraction x below 10^^. The effect 
of shadowing on the D meson Raa was estimated using the next-to-leading order (NLO) perturbative 
QCD calculation by Mangano, Nason, and Ridolfi (MNR) f46!| with CTEQ6M parton distribution func- 
tions HTI and the EPS09NLO parametrization B8l of their nuclear modification. The uncertainty band 
determined by the EPS09 uncertainties is shown in the left-hand panel of Fig. [H together with the av- 
erage D meson Raa- The shadowing-induced effect on the Raa is limited to ±15% for px> 6 GeV/c, 
suggesting that the strong suppression observed in the data is a final-state effect. 

The expected colour charge and parton mass dependences of parton energy loss should be addressed 
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Table 5: Average Raa as a function of pi for prompt D mesons in the 0-20% and 40-80% centrality classes. The 
systematic error does not include the normalization uncertainty, which is ±5.3% (±7.5%) for the 0-20% (40-80%) 
centrality class. 



Pt interval 
(GeV/c) 


Raa ± s 
0-20% centrality 


tat ± syst 
40-80% centrality 


2-3 
3^ 

4- 5 

5- 6 

6- 8 
8-12 

12-16 


0.51 ±0.10+°;^^ 
0.37 ±0.06to\l 
0.33±0.05+U;}? 
0.27 ±0.07+°;°^ 
0.28 ±0.04+°;°^ 
0.26 ±0.031°;°^ 
0.35 ±0.06^°:}° 


0.75±0.13+|Jj^ 
0.59±0.09+Jj;^^ 
0.55 ±0.07+°}^ 
0.54 ± 0.08+°: 
0.60±0.08+Jj;}^ 
0.66 ±0.08+°]^ 
0.64±0.16+|J:}^ 



Table 6: Average Raa as a function of centrality for prompt D mesons in the transverse momentum interval 
6<pt< 12GeV/c. 



Centrality 


/?AAi Stat ± syst(uncorr) ± syst(corr) 


0-10% 
10-20% 
20-10% 
40-60% 
60-80% 


0.23 ±0.03;°:°^!°:°^ 
0.28 ± 0.04 

42 + 04 +°°*^ 

0.81 ±0.10 



by comparing the nuclear modification factor of D and K mesons. Since final results on the pion Raa 
at the LHC are not yet available, we compare here to charged particles. Preliminary results 0491 have 
shown that the chai^ged-pion Raa coincides with that of charged particles above pt ~ 5 GeV/c and it 
is lower by 30% at 3 GeV/c. The comparison between D meson and charged particle Raa, reported in 
the right-hand panel of Fig. [H shows that the average D meson nuclear modification factor is close to 
that of charged particles |[T9l. However, considering that the systematic uncertainties of D mesons are 
not fully coiTclated with pt, there is an indication for 7?^^ > R^^^^^'^. In the same figure, the nucleai^ 
modification factor measured by the CMS Collaboration for non-prompt J/\j/ mesons (from B decays) 
with Pt > 6.5 GeV/c |[25l is also shown. Their suppression is clearly weaker than that of charged 
particles, while the comparison with D mesons is not conclusive and would require more differential and 
precise measurements of the transverse momentum dependence. 

Several theoretical models based on parton energy loss compute the charm nuclear modification factor: 
(I) MM, (11) USB, (III) El, (IV) (H, (V) El, (VI) ES, (VII) m, (VIII) dl. Figure i displays 
the comparison of these models to the average D meson Raa, for central Pb-Pb collisions (0-20%), 
along with the comparison to the charged-particle Raa El, for those models that also compute this 
observable: (I) US, (II) EU, (HI) [52], (VII) 1561. Among the models that compute both observables, 
radiative energy loss supplemented with in-medium D meson dissociation (I) |[T3l and radiative plus 
collisional energy loss in the WHDG (II) IBTl and CUJETl.O (VII) |[56l implementations describe rea- 
sonably well at the same time the charm and light-flavour suppression. While in the former calculation 
the medium density is tuned to describe the inclusive jet suppression at the LHC [50], for the latter two 
it is extrapolated to LHC conditions starting from the value that describes the pion suppression at RHIC 
energy (-y/%J^ = 200 GeV). This could explain why these two models are somewhat low with respect to 
the charged-particle Raa data. A model based on AdS/CFT drag coefficients (III) |[52ll underestimates 
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significantly the chami /?aa and has very limited predictive power for the light-flavour Raa- 
8 Summary 

The first ALICE results on the nuclear modification factor Raa for charm hadrons in Pb-Pb collisions at 
a centre-of-mass energy y/s^ = 2.76 TeV indicate strong in-medium energy loss for charm quarks. The 
D'', D+, and D*+ Raa, measured for the first time as a function of transverse momentum and centrality, 
is in the range 0.25-0.35 for 5 < pt < 16 GeV/c for the 20% most central collisions. For pt below 
5 GeV/c, and towards peripheral collisions, there is a tendency for an increase of Raa for D" mesons. 

The suppression is almost as large as that observed for charged particles, which are mainly light-flavour 
hadrons, with a possible indication, not fully significant with the present level of experimental uncer- 
tainties, of R^^ > R^^^^^'^. The expected effect of PDF nuclear shadowing is small (< 15%) above 
Pt = 6 GeV/c, indicating that the large measured suppression cannot be explained by initial-state effects. 
Some of the pQCD models based on various implementations of parton energy loss succeed reasonably 
well at describing simultaneously the suppression of light flavour and charm hadrons. 

The precision of the measurements will be improved in the future, using the large sample of Pb-Pb 
collisions recorded in 2011. In addition, p-Pb collision data expected in 2013 will provide insight on 
possible initial-state effects in the low-momentum region. 
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